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Abstract 

The clustering of ultra high energy cosmic ray events suggests that they 
have originated from compact sources. One of the possible physical mechanisms 
by which ultra high energy nuclei reach the Earth from far away astrophysi- 
cal sources (quasars or BL Lac objects) evading the Greisen-Zatsepin-Kuzmin 
(GZK) cut-off is by violation of Lorentz invariance. Assuming that there is 
violation of Lorentz invariance, we calculate the expected number of neutron 
events from some of the EGRET sources (including 7-ray loud BL Lac objects) 
which can be correlated in direction with ultra high energy cosmic ray events 
observed by AGASA above energy 4 x lO^^eV. We present in this paper what 
AGASA should see in future if violation of Lorentz invariance is responsible 
for the propagation of ultra high energy cosmic rays having energies above the 
GZK cut-off when there is a correlation of EGRET sources with the ultra high 
energy cosmic ray events. 
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Observational data on ultra high energy cosmic ray events have been recorded by 
different experiments like Fly's Eye [1, 2], Akeno Giant Air Shower Array (AGASA) 
[3], Yakutsk experiment [4], Haverah Park [5], Volcano Ranch experiment [6] and 
Sydney University Giant Air-shower Recorder (SUGAR) [7]. The composition of 
the ultra high energy cosmic rays (UHECR) has been studied in [8, 10]. Data from 
Fly's Eye experiment [8] suggest that the chemical composition is dominated by 
heavy nuclei up to ankle (lO^^'^)eV and thereafter by lighter component like protons. 
The AGASA data [10] suggest a mixed composition of both protons and heavier 
nuclei. The present experiments do not give us enough information about the chemical 
composition of primary UHECR. This topic has been addressed in [11]. If the UHECR 
are nuclei, then what could be their origin and how they propagate in the extragalactic 
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magnetic field are at present exciting fields of study. References [12, 13, 14] discuss 
different suggestions on these issues. 

The Greisen-Zatsepin-Kuzmin(GZK) cut-off [15] will prevent nucleons of energy 

greater than 4 x lO^^eV from travelling more than about 50 Mpc, but most of the 
suitable astrophysical acceleration sites are located at greater distances. AGASA ex- 
periment [3] has recorded about 8 events above 10^'^eV, whereas data recorded by 
HiRes fluorescence detector [16] seem to be in agreement with the GZK-cutofF [17]. 
Since currently the experimental data are not convincing enough for accepting the 
presence of the GZK cut-off, theoretical suggestions for observability of cosmic rays 
above the GZK cut-off remain of interest to us. The Physics involving the formation 
and annhilation of topological defects (TDs) has been studied as a possible explana- 
tion of observing UHECR. However TD annhilation has unique observational conse- 
quences such as, copious production of particles like neutrinos and 7-rays [18]. The 
experimental signatures of the Z-burst scenario has been discussed in [19]. Gamma 
ray bursts have also been identified as a possible explanation of observing UHECR 
above the GZK cut-off [20, 21]. 

If the highest energy cosmic rays have a neutral component then that will come to the 
Earth from the direction of its source. Elbert and Sommers [22] identified radio-loud 
quasar 3C147 as an ideal source within 10° of the highest energy event observed by 
Fly's Eye experiment [1] at energy 3.2 x 10^°eV. Later another quasar PG0117 -|- 213 
was identified by Biermann [23] within the error cone of the second highest energy 
event observed by AGASA [24]. In the last few years a lot of work has been done 
by different groups [25, 26, 27, 28] in correlating UHECR events with radio quasars 
and BL Lac objects. It was pointed out by Dubovsky et al. [29] that the statistics 
of clustering of ultra high energy cosmic ray events suggests their correlation with 
compact sources. In [30] the authors have pointed out that there is no significant 
correlation between quasars or BL Lac objects and ultra high energy cosmic ray 
events. They have used data from Haverah Park [5] and Volcano Ranch experiments 
[6] for their statistical analysise. However, the authors of [26] have shown that there 
is a significant correlation between compact sources and observed ultra high energy 
cosmic ray events using the data from AGASA [3] and Yakutsk [4] experiments. The 
probability of positional coincidence between BL Lac objects and UHECR events to 
occur by chance in a random distribution is of the order of 10^^. The still unknown 
data from HiRes experiment [2] and the accumulating data from AUGER experiment 
[31] will enable us to carry out statistical analysis with a bigger set of data in the 
near future. 

Since, as of now it is not clear whether the UHECR events are statistically uncor- 
related with compact sources at high redshifts, the physical mechanisms by which 



2 



neutral component of UHECR may come to the Earth evading the GZK cut-off from 
high redshifts are still exciting fields of study for us. In the Z-burst scenario [32] 
or in models of hypothetical "immune messengers" [33] neutral particles can come 
to the Earth from the direction of the source evading the GZK cut-off. It has been 
suggested by Coleman and Glashow [34] that if there is violation of Lorentz invari- 
ance (VLI) then protons and neutrons of energies above the GZK cut-off will reach 
the Earth without interacting with the cosmic microwave background. In the present 
work we show that one can test whether VLI is responsible for the observability of the 
UHECR events above the GZK cut-off using the data from AGASA experiment and 
certainly with the data from AUGER [31] experiment in future. We have predicted 
what AGASA should see in future if VLI is the underlying physical mechanism for 
propagation of cosmic rays with energies above the GZK cut-off from compact sources 
at high redshifts. 

We consider the ultra high energy cosmic ray events observed by AGASA above the 
GZK cut-off. Among this set of events we find that some of them can be correlated 
in direction with EGRET sources from the third EGRET catalog [36]. In [27] the 
authors have suggested that the sources of UHECR are high energy peaked BL Lac 
objects. Our set of correlated EGRET sources also include 7-ray loud BL Lac objects. 
If the radiation energy density is sufficiently high in a source, photo-pion production 
leads to the generation of a sufficient number of neutrons which can escape from the 
system. If there is VLI the neutrons having an energy above a certain energy will not 
decay. The energy above which the neutrons become stable depends on the degree of 
VLI. In [35] the authors have obtained a limit on the energy above which neutrons 
become stable, using observational data from the Yakutsk experiment [4]. These 
neutrons can travel through the cosmic microwave background with energy above the 
GZK cut-off when there is VLI and in that case we expect to detect them from the 
direction of the compact source. For clusters of UHECR events we consider sources 
within 4° of the events and for a single event we correlate with a source within 2.5°. 
In our correlated data set we have eight EGRET sources, among which three (3EG 
J0433+2908, 3EG J1052+5718 and 3EG J1424+3734) are confirmed BL Lac objects 
from the Veron 2001 catalog [37, 27]. Among the remaining five, 3EG J 1744-03 10 
is a quasar [38]. The four unidentified EGRET sources i.e. 3EG J1824+3441, 3EG 
J1903+0550, 3EG J0429+0337 and 3EG J0215+1123 may be BL Lac objects that 
have not yet been confirmed. Two of these four unidentified EGRET sources are also 
present in Table 3 of [27]. 
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TABLE I. The EGRET sources correlated with UHECR events above 4 x lO^'^ eV 
in AGASA data 



3EG J 


EGRET ID 


Type of 


Is 


hs 


z 


/ 


h 


E 






Object 


(deg) 


(deg) 




(deg) 


(deg) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


0433+2908 


AGN 


BL Lac 


170.5 


-12.6 


> 0.3 


170.4 


-11.2 


5.47 






2EG J0432+2910 


























171.1 


-10.8 


4.89 


1052+5718 


Possible AGN 


BL Lac 


149.6 


54.42 


0.144 


147.5 


56.2 


5.35 






RGB J1058+564 


























145.5 


55.1 


7.76 


1424+3734 




BL Lac 


63.95 


66.92 


0.564 


68.5 


69.1 


4.97 






TEX 1428+370 














1744-0310 


AGN 


Quasar 


22.19 


13.42 


1.054 


22.8 


15.7 


4.27 


0215+1123 






153.75 


-46.37 




152.9 


-43.9 


4.2 


0429+0337 






191.4 


-29.08 




191.3 


-26.5 


6.19 


1824+3441 






62.49 


20.14 




63.5 


19.4 


9.79 


1903+0550 






39.52 


-0.05 




39.9 


-2.1 


7.53 



In the first row of TABLE I we note that the redshift of the EGRET source 3EG 
J0433+2908 is more than 0.3 [39]. The fourth and fifth columns of TABLE I give 
longitudes and latitudes of the sources in Galactic coordinates and the seventh and 
eighth columns display longitudes and latitudes of the UHECR events in Galactic 
coordinates. 



We consider the sources of TABLE I whose redshifts are known and calculate 
the number of neutron events from them expected to be detected by AGASA in 30 
years. Wc calculate the expected number of neutron events from the source 3EG 
J0433+2908 assuming its redshift to be 0.3. In future we will come to know the 
redshifts of the other EGRET sources and then it will also be possible to calculate 
the expected number of neutron events from them. There are both theoretical [40] and 
observational [41] reasons to believe that when proton acceleration is being limited 
by energy losses, the luminosity of the object in very high energy cosmic rays Lqr 
is approximately equal to its luminosity in gamma rays. The cosmic ray luminosity 
LcR in the energy range 1 x 10^^eV< £■ < 4 x lO^'^eV can be assumed to be emitted 
equally in each decade of energy E. In that case, in each decade of energy we expect 
the power emitted in UHECR in the energy range of 1 x lO^^eV to 4 x 10^°eV to be 
approximately L^/j/lO. Let 'A' be the area of AGASA detector which is IQ^'^cm? [3]. 
Here we mention that the exposure of the AGASA detector is energy independent in 
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the energy range in which we are interested. One can see the plot of exposure against 
energy of the UHECRs for AGAS A detector in [42] . The expected number of neutron 
events in a time interval dto in AGASA within the source energy interval of Ei and 
E2 can be expressed as 



where d is the luminosity distance of the source from the Earth. N'^ and are 
respectively the number of neutrons emitted at the source and observed by the detec- 
tor. Es, Eo are the source energy and observed energy of the neutrons respectively. 
If the redshift of the source is z then E^ — Es/{1 + z). Similary the correction to be 
applied to observed time to due to redshift of the source is tg — ts{l + z), where tg is 
the time at which the neutron is emitted from the source. 

We have assumed there is VLI and that the lower limit of the above integration 
is such that neutrons are stable above this energy. 

We define F^'' as the energy received in gamma rays per second per cm^ on the 
surface of the Earth from an EGRET source and this quantity can be calculated from 
[36]. e„ is the efficiency of neutron production in the source. When the photon density 
in a source is sufficiently high then the efficiency of neutron and proton production 
become comparable near the end of the UHECR spectrum [43]. Hence it is not 
unreasonable to assume e„ = 1/2 in our calculation. If we assume the luminosity of 
the source in gamma rays of energy more than lOOMeV L^, to be a fraction x times 
the luminosity of the source in UHEGR Lcr, then eqn.(l) can be written as 



An almost similar procedure has been followed in ref.[43] in calculating number of 
neutron events from Centaurus A but without the assumption of VLI. Using eqn. (2) 
we calculate the expected number of neutron events between observed energy intervals 
of 1 X lO^^eV and 2 x lO^^cV in AGASA from the sources of TABLE 1 whose redshifts 
are known. The second column of TABLE II displays the photon flux in 10^^ photon 
cm~'^sec~^ above 100 MeV energy for the EGRET sources in the first column [36]. 
The third column shows the photon spectral indices from [36] of the EGRET sources 
displayed in the first column. The fourth and fifth columns present the expected 
number of neutron events in a time interval of 30 years in AGASA for case (1) 
e„ = 1/2, X = 1 and for case(2) e„ = 1/2, x = 1/2 respectively. 




(1) 




(2) 
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TABLE II. The number of neutron events between energy 1 x 10^° eV and 2 x lO^^eV 
expected in AG AS A in 30 years from some of the EGRET sources of TABLE I. 



3EG J 


Photon flux 


Photon Spectral 


Number of Neutron 


Number of Neutron 




above 100 MeV 


Index 


Events in case (1) 


Events in case (2) 


0433+2908 


22.0 


1.9 


8.4 


4.2 


1052+5718 


6.5 


2.51 


2.2 


1.1 


1424+3734 


16.3 


3.25 


3.5 


1.7 


1744-0310 


21.9 


2.42 


4.29 


2.1 



In the last 10 years data from AGASA there is no event above lO^^'eV energy from 
the direction of the sources of first column of TABLE II. If VLI is the underlying 
mechanism for UHECR propagation then we can expect to see UHECR events above 
10^°eV from the direction of the sources presented in TABLE II in future. If we 
increase the area of the detector or the time of data collection then of course we 
would detect more neutron events than that presented in TABLE II. The number of 
neutron events will increase linearly with increasing area or time of data collection 
by the detector. Since the area of AUGER is 30 times larger than AGASA, we can 
expect such neutron events in a much shorter interval of time in AUGER. 



Conclusion 

In this paper we have presented how the UHECR data above the GZK cut-off can 
tell us whether VLI is the underlying mechanism for the propagation of these cosmic 
rays from quasars or BL Lac objects at high redshifts. 



6 



References 

[1] D. J. Bird et al., Astrophys. J. 441, 144 (1995). 

[2] T. Abu-Zayyad et al.[High Resolution Fly's Eye Collaboration], arXiv:astro- 
ph/0208243. 

[3] M.Takeda et al, Phys. Rev. Lett. 81, 1163 (1998); M.Takeda et al., arXivrastro- 
ph/9902239; N. Hayashida et al., Astrophys. J. 522, 225 (1999); M. Takeda 
et al., Astropart. Phys. 19, 447 (2003); see also http://www-akeno.icrr.u- 
tokyo. ac.jp/AGASA/. 

[4] B. N. Afanasiev ct al., Proc. Int. Symp. on Extremely High Energy Cosmic Rays: 
Astrophysics and Future Observatories (ed. M. Nagano, Institute for Cosmic Ray 
Research, Univ. of Tokyo) p. 32 (1996). 

[5] M. Ave, J. A. Hinton, R. A. Vazquez, A. A. Watson and E. Zas, Phys. Rev. Lett. 
85, 2244 (2000): M. Ave, J. Knapp, J. Lloyd-Evans, M. Marchesini and A. A. 
Watson, Astropart. Phys. 19, 47 (2003). 

[6] J. Linsley, Catalog of Highest Energy Cosmic Rays, (World Data Center of Cos- 
mic Rays, Institute of Physical and Chemical Research, Itabashi, Tokyo)p.3 
(1980). 

[7] M. M. Winn, J. Ulrichs, L. S. Peak, C. B. Mccusker and L. Horton, J. Phys. G 
12, 653 (1986). 

[8] D. J. Bird et al, Phys. Rev. Lett. 71, 3401 (1993); Astrophys. J.424, 491 (1994). 
[9] T. Abu-Zayyad et al., astro-ph/9911144; Astrophys. J. 557, 686 (2001). 
[10] N. Hayashida et al., J. Phys. G 21, 1101 (1995). 

[11] B. R. Dawson, R. Meyhandan and K.M.Simpson, Astropart. Phys.9, 331 (1998). 
[12] F. W. Stecker, Phys. Rev. Lett. 80, 1816 (1998). 

[13] Eun-Joo Ahn, G. Medina- Tanco, P. L. Biermann and T. Stanev, arXiv: astro- 
ph/9911123. 

[14] G. Bertone, C. Isola, M. Lemoine, G. Sigl, Phys. Rev. D 66, 103003 (2002). 

[15] K. Greisen, Phys. Rev. Lett. 16, 748 (1966); G. T. Zatsepin and V. A. Kuzmin, 
Sov. Phys. JETP Lett. 4, 78 (1966). 



7 



[16] C. H. Jui (HiRes collaboration), Proc. 27th International Cosmic Ray Confer- 
ence, Hamburg 2001. 

[17] J. N. Bahcall and E. Waxman, Phys. Lett. B 556, 1 (2003). 

[18] G. Sigl, Space Sci. Rev. 75, 375 (1996) and refs. therein. 

[19] T. J. Weiler, Astropart. Phys. 11, 303 (1999). 

[20] E. Waxman, arXiv:astro-ph/0210638. 

[21] M. Vietri, D. D. Marco and D. Guetta, Astrophys. J. 592, 378 (2003). 

[22] J. Elbert and P. Sommers, Astrophys. J. 441, 151 (1995). 

[23] P. Biermann, Nucl. Phys. (Proc. Suppl.) B 43, 221 (1995). 

[24] N. Hayashida et al., Phys. Rev. Lett. 73, 3491 (1994). 

[25] G. R. Farrar and P. L. Biermann, Phys. Rev. Lett. 81, 3597 (1998). 

[26] P. G. Tinyakov and I. I. Tkachev, JETP Lett. 74, 1 (2001); JETP Lett. 74, 445 
(2001); Astropart. Phys. 18, 165 (2002); 

[27] D. S. Gorbunov, P. G. Tinyakov, I. I. Tkachev and S. V. Troitsky, Astrophys. J. 
577, L93 (2002). 

[28] A. Virmani, S. Bhattacharya, P. Jain, S. Razzaque, J. P. Ralston and D. W. 
McKay, Astropart. Phys. 17, 489 (2002). 

[29] S. L. Dubovsky, P. G. Tinyakov and I. I. Tkachev, Phys. Rev. Lett. 85, 
1154(2000). 

[30] D. F. Torres, S. Reucroft, O. Reimer and L. A. Anchordoqui, Astrophys. J. 595, 
L13 (2003). 

[31] http://www.auger.org/ 

[32] D. Fargion, B. Mele and A. Salis, Astrophys. J. 517, 725 (1999). 

[33] D. J. Chung, G. R. Farrar and E. W. Kolb, Phys. Rev. D 57, 4606 (1998). 

[34] S. Coleman and S. L. Glashow, Phys. Rev. D 59, 116008 (1999) ; arXiv: hep- 
ph/9808446. 

[35] S. L. Dubovsky and P. G. Tinyakov, Astropart. Phys. 18, 89 (2002). 

8 



[36] R. C. Hartman et al., Astrophys. J. Supp. 123, 79 (1999). 

[37] M. -P. Veron-Cetty and P. Vcron, Astron. Astroph. 374, 92 (2001). 

[38] A. Hewitta and G. Burbidge, Astrophys. J. Supp. 87, 451 (1993). 

[39] J. P. Halpern, M. Eracleous and J. R. Mattox, Astron. J. 125, 572 (2003). 

[40] K. Mannheim, Astron. Astrophys. 269, 67 (1993). 

[41] S. Rawhngs and R. Saunders, Nature 349, 138 (1991). 

[42] M. Takeda et al, Astropart. Phys. 19, 447 (2003). 

[43] L. A. Anchordoqui, H. Goldberg and T. J. Weiler, Phys. Rev. Lett. 87, 081101 
(2001). 



9 



